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Evaluation  of  Singular  Electric  Field  Integral  Equation  (EFBE)  matrix  elements 
Introduction 

The  electric  field  integral  equation  (EFIE)  solution  for  scattering  from  an  arbitrary 
three-dimensional  geometry  was  outlined  in  a  paper  by  Rao,  Wilton  and  Glisson  (RWG) 
in  1982  [1],  Since  that  time  their  formulation  has  been  incorporated  into  many  computer 
codes  to  solve  a  large  variety  of  scattering  and  radiation  problems,  e.g.  EIGER, 
CarlosSD,  FISC,  Feko,  MoM3D,  IBC3D  [2-7].  Experience  has  taught  us  that  the 
accuracy  of  quantities  derived  from  the  original  RWG  solution  are  a  function  of  the 
electrical  size  of  the  problem,  which  is  reflected  in  the  order  of  the  system  of  equations 
that  must  be  solved,  and  the  electrical  dimensions  of  the  triangular  patches  [8], 

The  ill-conditioned  matrices  that  result  from  electrically  small  patches  can  be 
remedied  using  of  the  ‘loop-star’  approach  [9].  For  electrically  large  problems,  with 
large  order  systems  of  equations,  the  physical  fidelity  of  the  solution  depends  on  the 
numerical  accuracy  of  the  Galerkin  inner  products  which  populate  the  method  of 
moments  (MoM)  impedance  matrix.  A  new  treatment  of  these  inner  products  is 
presented  here  which  improves  the  accuracy  of  the  singular  surface  integrals  that  result 
from  the  Galerkin  formulation. 

The  question  of  accuracy  has  been  a  subject  of  research  in  the  computational 
mechanics  field  for  some  time.  Recently,  the  computational  mechanics  community  has 
begun  formulating  boundary  integral  equation  (BEE)  solutions  of  Fredholm  integral 
equations  using  a  Galerkin  descretization  [10,11],  The  focus  of  their  research  was  to 
obtain  high  accuracy  in  the  numerical  treatment  of  the  resulting  singular  4-dimensional 
integrals.  The  methods  developed  for  computational  mechanics  can  be  directly  applied  to 
the  Galerkin  (RWG)  solution  of  the  EFIE  to  obtain  high  accuracy  in  the  numerical 
integration  of  the  singular  surface  integrals  that  are  fundamental  to  the  numerical 
solution. 

Galerkin  Solution  of  Electromagnetic  EFIE 

The  traditional  RWG  solution  will  be  summarized  here  for  objects  that  are 
represented  by  perfectly  conducting  surfaces.  The  time  dependence  is  given  by 
The  electromagnetic  problem  is  formulated  using  the  electric  field  integral  equation 

(EFEE)  which  is  then  used  to  compute  unknown  electric  currents,  J ,  on  the  surface  of  the 
object.  These  surface  currents  are  integrated  to  produce  the  desired  quantity,  for 
example,  the  far  field  scattering  amplitude. 

The  unknown  currents  are  expanded  in  terms  of  RWG  vector  basis  functions  on 
the  surface.  A  single  RWG  basis  set  element,  shown  in  Figure  1,  consists  of  a  pair  of 
triangular  patches  that  share  a  common  edge.  The  Galerkin  solution  of  the  EFIE 
generates  a  system  of  linear  equations  which  has  the  amplitudes  of  the  RWG  current 
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expansion  as  unknowns.  This  system  of  equations  can  be  written  in  the  following 
compact  form, 


ZJ  =  F.  (1) 

The  square  matrix  Z  is  dense  and  complex  valued.  This  system  of  equations  is 

solved  using  standard  numerical  linear  algebra  methods.  The  vector  V  represents  the 
electric  field  excitation  on  the  structure,  usually  a  plane  wave.  The  individual  elements 
of  the  Z  matrix  are  inner  products  of  an  RWG  basis  and  testing  function  vdth  the  EFIE 
operator  for  the  PEC  surface  boundary  condition.  These  inner  products  generate  double 
integrals  over  pairs  of  RWG  basis  functions  and  appear  as  a  sum  of  4  terms, 

(2) 

p,q^±l 

Each  term  in  the  sum  represents  the  interaction  of  one  triangle  with  another  in  the  RWG 
basis  and  testing  function  pair.  This  interaction  is  computed  using  the  following 
expression. 


fpypi 


1 


g(F,F'), 


where  g{r,r')  is  the  free  space  Greens  function. 


gir,r') 


4;r|r-r' 


(3) 


(4) 


and. 


,1^0  -  Speed  of  light,  /  =  wave  frequency. 


1]  =  ,1— ,  impedance  of  free  space, 

Aj,Al  =  area  of  each  triangle, 
p%pt  =  vector  from  vertex  opposite  edge  j,k, 
Sj,S^  =  ±1,  current  flow  toward  or  away  from  edge, 
/.,4  =  edge  lengths  associated  with  RWG  triangle  basis. 


A  numerical  difficulty  arises  in  evaluating  (3)  over  regions  where  the  denominator  of  the 
Greens  function  becomes  zero.  When  this  occurs,  the  weak  singularity  in  the  integrand 
of  the  integral  must  be  interpreted  in  a  Cauchy  principle  value  sense. 
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The  treatment  of  this  singularity  has  generally  consisted  of  employing  a 
singularity  extraction  approach  where  the  Greens  fxmction  in  the  integrand  is  split  into 
two  parts. 


g{7,r) 


4;rlF-r' 


4;r|r  -  r'  4;r|F  -  r ' 


(5) 


This  produces  two  integrals,  the  first  has  a  regular  integrand,  which  is  numerically 
integrated,  and  the  second  that  is  evaluated  analytically  [12],  This  approach  is  limited 
because  the  integrand  still  contains  a  R~'^  term  to  be  integrated  numerically,  requiring  a 
large  number  of  evaluations  for  a  given  accuracy.  This  method  has  also  been  applied  to 
formulations  with  curved  geometry  where  the  singularity  extraction  is  the  first  term  in  a 
Taylor’s  series,  however  this  approximate  method  can  be  improved  upon  with  the 
formulation  outlined  here. 

Erichsen  and  Andra  [10,11]  encountered  integrals  of  the  type  shown  in  (3),  and 
others  with  higher  order  singularities,  and  were  able  to  completely  remove  the  singularity 
using  a  combination  of  relative  coordinates,  changing  the  order  of  integration  and 
Duffy’s  transformations  [13].  In  addition  to  completely  removing  the  singularity  their 
method  also  allowed  some  parts  of  the  4-dimensional  Galerkin  integrals  to  be  evaluated 
analytically,  thereby  reducing  the  effort  required  to  numerically  evaluate  the  entire  inner 
product.  Their  technique  for  regularizing  multi-dimensional  singular  surface  integrals  is 
applied  here  to  the  traditional  RWG  EFIE  solution. 


Regularization  of  Singular  Surface  Integrals 

The  integrals  in  (3)  can  be  categorized  into  5  classes.  The  criterion  is  based  upon 
the  geometrical  relationship  between  the  two  triangular  regions,  labeled  Af,A^.  The 

fundamental  distinction  is  whether  or  not  the  two  triangles  intersect  in  some  manner,  and 
if  they  do  not,  then  how  close  or  far  apart  are  the  two  triangles.  The  5  classes  are 
summarized  in  the  following  table. 


Type 

Relationship 

Distribution  of  Singularity 

1 

Common  facet 

two  dimensional  planar 

2 

Common  edge 

one  dimensional  line 

3 

Common  vertex 

single  point 

4 

Distance  <  Facet  Dimension 

Near  Singular 

5 

Distance  >  Facet  Dimension 

Non- Singular 

Types  1  through  3  have  singularities  in  the  integrands  and  explicit  formulas  will 
be  developed  to  analytically  remove  the  singularity.  Type  4  integrals  are  described  as 
near  singular  and  are  studied  elsewhere  [14].  Type  5  present  no  singular  behavior  and 
can  be  evaluated  using  standard  numerical  methods  for  integrating  regular  functions  over 
a  triangular  region  [15]. 

The  procedure  for  transforming  (3)  into  a  regular,  non-singular,  surface  integral 
involves  the  following  steps. 
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1.  Introduce  relative  coordinates. 

2.  Represent  the  domain  of  integration  such  that  the  outer  integration  is 
over  the  relative  coordinates. 

3.  Recombine  symmetric  integration  domains. 

4.  Use  Duffy’s  coordinates  to  remove  the  weak  singularity  in  the 
remaining  singular  integrals. 

5.  Employ  ordinary  gaussian  quadrature  on  the  remaining  regular 
integrals. 

Simplex  Coordinates 

The  integral  in  (3)  is  evaluated  by  transforming  from  ordinary  3-dimensional 
Cartesian  coordinates  to  a  simplex  coordinate  system.  In  simplex  coordinates  each 
triangle  is  mapped  to  a  standard  unit  triangle.  The  unit  triangle  is  defined  by, 

(6) 

A  triangle  in  3-dimensional  Cartesian  space  with  vertices,  Vi,V2,V^,  is  mapped  onto  this 
unit  triangle  such  that, 


F,=^(^„4)  =  (0,0) 

V2^i4a2)=m  (7) 


Any  point  in  the  triangle  defined  with  three  vertices  can  now  be  described  using  the  2 
simplex  coordinates. 

Ignoring  the  constants  outside  the  integral  shown  in  (3)  the  transformed  integral 
using  two  sets  of  simplex  coordinates  is, 


1 


1  vi 


g(f,r)dA;dA1^2ApAl\  I  I  }  (8) 


^1  =0  ^2  ^  ^  ^ 


where  the  integrand  has  now  been  written  as  ratio  of  a  regular  function  divided  by  a 
distance  function.  The  constants  in  front  of  the  new  integral  are  the  Jacobian  of  the 
transform  from  Cartesian  to  simplex  coordinates. 


Decomposition  of  the  domain 

The  4-dimensional  Galerkin  inner  product  integral  is  decomposed  into  six 
separate  4-dimensional  integrals  in  the  process  of  converting  to  relative  coordinates  and 
interchanging  the  order  of  integration.  Explicit  details  of  the  process  are  presented  here 
because  the  method  allows  for  variations,  for  example  curved  triangular  patches  or 
rectangular  patches,  to  be  analyzed  using  the  same  approach. 
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The  fundamental  quantity  to  evaluate  is  the  Galerkin  inner  product  integral. 
Using  simplex  coordinates,  ,  the  integral  is  converted  to  the  follotving  standard 
form. 


j  f  }  ]  (9) 

5=0  ^j=Ov,=o  1/2=0  ) 

The  order  of  integration  for  ^2  and  t^i  is  interchanged  in  preparation  of  a  conversion  to 
the  relative  coordinates. 


1  1 


fill 

=0  7/1  =0^2=0  72  ==0 


/(£,7) 


(10) 


Relative  coordinates  are  introduced  next,  the  relative  coordinates  are  formed  from  the 
simplex  coordinates  in  the  following  manner. 


^2  ~  Vl  ^^2- 


(11) 


The  integral  using  relative  coordinates  as  the  integration  variables  is  now  written  in  the 
following  form. 


1  1-5  5  "1+5-52 

JIM 


^l=0«i^-fl^2=0  «2=-^2 


(12) 


The  domain  of  integration  of  this  integral  is  shown  in  Figure  2.  The  goal  is  to 
interchange  the  order  of  integration  in  the  integral  so  that  the  outer  integrations  are  over 
the  relative  coordinates  u\  and  ui. 

Examining  the  integration  domain  in  Figure  2  facilitates  the  interchange  of  order 
of  integration.  By  inspection  the  ui  integral  upon  interchange  becomes, 

I  j  \-di,du,+  \''l-di,du,.  (13) 

Note  that  the  original  domain,  a  single  trapezoid,  has  now  been  divided  into  two 
triangular  regions.  The  ^2,  U2  integration  interchange  produces  three  new  integrals, 

5  “1+5-52  0  5  «i  u,+^i  5+«i-<i2 

I  du^d^^-  I  J  d^^du^-v  I  ^■■•d^2^U2.  (14) 

52=0  «2  =-52  «2=-5  52=-«2  112=0^2=0  «2=“l  52=0 


The  original  domain  of  integration  has  been  divided  into  two  triangular  and  one 
rectangular  region. 
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Collecting  together  the  six  new  integrals  that  result  from  the  initial  interchange  of 
integration  variables  produces  the  following  set  of  terms, 


q  i, 


I  J  1  i 


(15) 


*'=  1  1  n 

U,=-l  f,=-U|  1(1=0  #2=0 


(16) 


0  1  U|+#,  #i-(i*2-“i) /■/ r  T7\ 

■*=  n  I  1 


#2=0 


«(#.(() 


(17) 


1  1-«|  0  #1  /'/  t  77\ 

=  i  J  J  1 

U,=0#,=0«2=-#i#2=-«2 


(18) 


1  l-U,  <1,  #1 


l/l=0  #1=0«2=0  #2=0  ■ 


(19) 


1  l-«,  «,+#! 


1  I  T  f  47|^«2''VfA. 

Ii,=0#i=0tt2=Wl  #2=0 


(20) 


The  next  step  is  to  interchange  the  order  of  the,  ui,  integration  in  the  six 
integrals.  Starting  with  l\  the  interchange  results  in  two  terms, 


1  0 


"l  1 


0  1 


f  J  l-di,du,+  J 

#I=-I/|  H2=-#,  K2=-1#,=-U2  1(2=1/,  #|=-U, 

The  completed  interchange  of  integration  variables  now  produces  a  split  l'  term, 

I'  =  ll+Il=l  I  1  1 


(21) 


«(#,(/) 

+  j  1  J  J 


IIj=-U2=-1#=-M2#2  =  “«2 

0  0  1#, 


(22) 


The  result  for  is, 
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‘^=1  n 

U,=-1«J=0^,=-H,  fj=0  ’“) 


is  similar  to  l’  because  it  produces  two  terms. 


0  0  1  -f  (  K  — \ 

i>  =  i|+r>  =  I  }  j  I 

to=-l  £,=0 


M2=Mi  ^2=0 

0  l+«i  1 


u  1+Hi  1  ginw2“«i^  r(  y.  yrx 

+  J  J  j  I 

u,=-lu,=0£=tt— u.  £.=0 


^=-l“2=Ofi=“2-“l  #2=0 


Interchanging  the  order  of  integration  for  the  remaining,  f*,  and  integrals  produce 
only  a  single  term  in  each  case. 


1  0  l-«,  f, 

•‘=1 J  1  1 

ttl=0  tt2=tti-l  #j=“«2  ^2=“‘«2 


/({.gx 

7J«.g) 


1  Wi  1-Wi 

=  1  n  J  . 

B|=0«2=0#,=0f2'=0 


1  0  l-a,  ^,-(u2-«,) 


«l=0  tt2=Wi  42"^ 


The  final  step  is  to  combine  integrals  together  with  overlapping  domains.  The  naming 
convention  here  is  used  to  reproduce  the  work  of  Andra, 


0014  0  #l-(«2-"l) 

1 1  1 U  +  ^  I 


«, — 1  -«i  4i — "1  t  ^2  ""“2  42  -4i  42 


R(4,u) 


i.i  1  J 


E,  =  .-j  J71 

wi=0 1/2=0  ^1=0  ^2=0 


0  I+«l  1  ^l-(«2-«l) 


E>  =  I?=  J  I  J  \ 

U,=-lM,=0  £=«.-!/,  f.=0 


l/l=“l  M2  =0  4\  =^2  -Wj  ^2  =0 


1  0  1-M, 


E.  =  r  =  J  J  j  J 

«J  =:0  Mj  =«!  -1  4^1  =-  ttj  ^2  =“«2  ^  ^  ^ 
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0  u.  I 

=  ‘!=f  1  I  J 


1  I 

E«  =  '‘=J  1  i  I 

ttl=Ou2=«l^I=W2-«l  ^2^ 


(32) 

R{^,u) 

(33) 

These  six  integrals,  Ei  through  E^,  represent  the  domain  decomposition  of  the  original 
single  domain  integral. 


Common  Facet 

The  Galerkin  inner  product  for  this  type  involves  integrating  over  identical 
triangular  facets.  The  evaluation  of  the  inner  product  integral  depends  on  the  functional 
form  of  the  singularity  in  simplex  coordinates.  The  Greens  function  spatial  distance 
term. 


=  (34) 

expressed  in  the  transform  coordinates  becomes, 


r=(l-^0V.+(^i-4)V2+^,V3 

7 -r  =  {11,-  )  V,  )%  +  (^,  -  %  (35) 

=  u,%-\-{u^-u,)%-uy^ 

and  is  a  function  of  the  relative  coordinates  only.  The  scalar  distance  is, 

R  =  Vk(V,  -  V,  -  V,)  +  u,%)-  (u,(V,  -  V,  -  V,)  +  u,% ), 

- -  (36) 

=  ■\ja,uf  +  2a2U,U2  +  a^ul , 

and  can  be  expressed  in  terms  of  a  simple  polynomial  in  the  relative  coordinates  only. 
The  coefficients  are. 


«>=(V,-V,-V3)-(V,-V2-V3), 

«2=(V,-V3-V3)-V3,  (37) 


This  allows  the  Greens  function  portion  to  be  moved  outside  the  inner  two  integrals  in  the 
inner  product. 
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The  remaining  parts  of  the  integrand  involve  dot  products  of  vectors  in  the  plane 
of  the  triangle,  each  vector  is  defined  with  a  starting  point  in  a  prescribed  vertex  of  the 
triangle  and  an  ending  point  that  resides  in  the  triangle.  Even  though  the  case  considered 
here  is  for  identical  facets,  or  triangles,  the  two  vectors  in  the  dot  product  are  not 
necessarily  the  same.  In  3-dimensional  Cartesian  space  let  these  two  vectors  be  denoted 

by. 


Pj=7-vj, 

A 


(38) 


Expressing  these  vectors  in  terms  of  the  simplex  coordinates  for  the  most  general  case, 
p,  =  r  -  V,  =  (1  -  -  ^2)%  +  V3  -  V,, 

A  =  =  (l-/7i)V;+(7i  -%)Vi  +T1X-V,. 


The  dot  product  generates  the  following  expression, 

pj-h  =  ((i-f,)v,  +«  -^2)v,+f,Vj-v,H(i-'/,)v;+(7,  -%)vi  +i7,v;-v.)  (40) 

Expanding  this  and  simplifying  the  resulting  expression, 

p,-p„=  ((l-^,)V,  +(#,  +f,V,-v,>(a-.;,)V,'+(>?,  -%)Vi  +7,V3'-v.) 

=  fe(v,+v,-v,)-f,v,+v,-vj-fe,(v'+v'-v,')-i7,v;+v;-v.)  (4i) 

~  P’i^inx  Pa^2^2  AA  PlVl  A^2  A 

where  again  the  function  can  be  written  as  a  polynomial.  The  coefficients  of  the 
polynomial  are  given  by, 

A-(v.+V3-Vi)-(V'+v'-v;) 

A  =  -(V2+V3-v,)-v' 

A  =  (V2+V3-v,)-(v;-fy) 

Pa  ”  p2 

A  =  V3-V'  (42) 

A  =  -V2-(v;-fy) 

A  =  (v,-v,)-(v'+v^v;) 

A  =  -(Vi-v,)-v' 

A  =  (v,-v,)-(v;-fy) 


The  same  dot  product  in  relative  coordinates  is. 


Pj-Pk=  +^,)+M("!+ft)+ A4(«. +f,)+ A("2  +&)fe 

+AS  +  Afe  +A(«,  +^,)+ A("2+A)+ A 

Collecting  together  the  expressions  for  the  Greens  function  and  the  dot  product  the 
original  integrand  becomes, 

^^4^*  S{rj')  =  ^h{^,u)g{u).  (44) 


It  is  this  function  that  will  be  used  as  the  integrand  for  the  six  integrals  that  span 
the  domain  of  integration  of  the  original  inner  product.  However,  the  six  integrals  can  be 
paired  together  to  form  three  integrals  because  of  the  symmetry  of  the  integration 
domains.  This  pairing  of  integrals  with  symmetric  domains  will  simplify  the  problem. 
First  consider  the  integrals  Ei  and  E2, 

0  0  1  ^i-(«2~«l)  1  _ 

^■=111  I  jK?,u)gm4:d4,dy^u„  (45) 

«,=-l«2=Wi 

^2=11  f  1  \>>(4,u)g{uyi4,d4,du^u,.  (46) 

«i  =0 1/2  =0  f^i  =0  ^2  =0 

Let  ^  =^'-u  and  u  =  -z  in  Ei.  The  results  are, 

0  0  1  ^  _ 

J  I  J  J 

0  0  1  ^j4-(Z2-2|)  .| 

=  111  1  (47) 

2l=l22=Zi#j=Zi  ^2=^2 

=  1  I  i"  1  ~K4' 

2,=0  22=0^,'=0^^=0^ 

The  limits  of  integration  are  identical  to  that  of  E2  so  that  the  two  integrals  can  be 
combined  together, 

Ji=E, +£3=1  f  J  '^}-h(^,u)g{u)  +  Jh{^  +  u-u)g{-u)d^^d^^du^du^.  (48) 

U|=OU2=0^,=0^2=o'^  ^ 

The  Greens  function  is  an  even  function  of  the  argument, 

g(-w)  =  g(w),  (49) 


the  new  integral,  J\,  becomes. 
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1  a,  1-a,  #i  I  _  _ 

J,=E, +£2=  III  ^  ,u)  +  h{^  -{■u-u))^{u)d^2'‘^^^du2du■^.  (50) 

Bi=0Uj  =0^1  =0^2=0^ 

The  integrals  E3  and  E4  can  be  combined  together  in  a  similar  fashion.  Substitute 
^  ^ «  and  u  =-z  in  E3, 

0  1+Ui  I  1 

J  J  I  1 

“l=-I  "2=0  ^l=“2-“l  6=0  j  ^ 

J  J  11  + 

2j=0  22“2i“1 

This  is  combined  with  E4, 

=  1  J  f  1  (52) 

|/j=K)tt2'«i-l  ^1=“M2  ^2=“«2 


to  produce, 

J,=E3  +  E,=  |  j  I  ]  ^^(^.u)  +  h(i+u-u)}g(u)d4^4,du,du,.  (53) 

«j=0«2=«i-l  ^i=-M2  ^2=-“2 


Lastly,  combine  E5  and  Ee  together  by  substituting  ^  =^'-u  and  u  =-z  in  E5, 


1  1  l-Oi  f|-(a2-«i)  1  ,  _  _ 

J3  =  E3  +  E,=  JJ  I  j  +  (54) 


l/j  =0  «2  ^2 


At  this  point  the  original  integral  has  been  successfully  transformed  into  a  sum  of  three 
integrals  with  the  desired  order  of  integration,  this  result  can  be  written  simply  as. 


1  6  I 


f  }  f 

6=06=0;?i=0,2=0^(^’'7) 


(55) 


;=1 


However,  the  singularity  in  the  Greens  function  is  still  present.  The  next  step  is  to 
employ  a  series  of  Duffy  coordinate  transformations  that  will  eliminate  the  singularity. 
The  removal  of  the  singularity  will  also  convert  the  three  integrals  into  one  regular 
integral. 

Duffy  Transformations 

The  Duffy  coordinate  transform  is  used  to  remove  a  singularity  at  the  origin  in 
two  dimensional  and  higher  order  integrals.  Each  of  the  J  integrals  has  a  singularity  at 
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the  origin,  namely  w  =  0 .  The  (wi,W2)  coordinate  integration  domains  as  shown  in  Figure 
3  for  each  of  the  Ji  integrals. 

In  each  of  the  J  integrals  the  Greens  function  contains  a  singularity  that  is  a 
function  of  the  relative  coordinates  only.  In  relative  coordinates,  the  Greens  function  for 
the  common  facet  configuration  is. 


Jkyja,Ui  +2a2U|M2  +a'3«2 

g(u)  = - ,  (56) 

+2«2W,M2 +  ^3^2 

The  Ji  integral  over  the  relative  coordinates  uses  the  standard  two-dimensional  Duffy 
transform.  Consider  the  part  of  Ji  that  is  the  integral  of  the  Greens  function  over  the 
relative  coordinates  only, 

1  »i  Jk^aiuj+la^UiU-^+a^ul 

J  f - . —  -.^"•du^du^.  (57) 

„,=o«;=o4;Ti/a,wf  -h2«2WiM2  +cCiUl 


Make  the  substitution. 


(58) 


of  these  coordinates  into  the  integral, 

I  “l  rir  +2o2«i«j  +0^1  •  I  jk(o-Jayi^a^+a^ 

II  .  .  r  r  . 


n  JkyjaiU^-i-2a2UiU2+a^U2  1  1  U'kco 

- p-  ^  -du^du,  =  f  J - ^ 

U|=0  ii2=0  42r-Y  OTjZ/i  “f”  2o^2^|1^2  ^2^2  oj=0rj=0^^^^^\ 


+  2^27  + 


(O'-’dridco . 


(59) 

The  singularity  has  been  converted  into  a  single  term,  <y,  in  the  denominator  that  is 
exactly  canceled  out  by  the  Jacobian  of  the  coordinate  transformation  in  the  numerator. 
The  resulting  non-singular  integral  is. 


1  I  jkco\lai+2ajTj+a^T7^ 

J,  =  n  ^  -drjdo,.  (60) 

«=0  n=o  4^  V  a^  +  2«27  +  O',;; 

However,  J2  is  not  in  the  proper  form  for  a  direct  Duffy  coordinate  transform, 

1  0 

J2  ~  J  J’**  5 

«[  =0  I/2=Mi-1 

but  can  be  put  converted  to  standard  form  by  substituting  2  =  w,  -  w, , 


(61) 


Evaluation  of  Singular  Electric  Field  Integral  Equation  matrix  elements 


13 


i  1 

J2=  I  \--dzdu^. 

aj=0  z=Mi 


The  Duffy  transform  for  this  type  is, 


W,  =  7J(0, 

z  =  0), 


(62) 


(63) 


converting  back  to  the  original  relative  coordinates  produces  the  necessary  transform. 


U^  =  T}0), 
Z-U^-U2=CO, 

Wj  =«j-6)  =  <y(;7-l). 


(64) 


Using  this  transform  in  the  integral  Ja, 

JkT]a,ul  +2a2"i<'2 


1  1 

11 

17=0 


f  f  - -  "-du^du^  = 

«r=o»j=Vi  47^^Ja^uf  +  2a2U^U2  + 

)+a3<»^  (7-1)^ 

4;r^«j  +  2a2}j(0^(jj  - 1)  +  a^co^(7j  - 1)^ 


(65) 


0)---d7]da>, 


removes  the  singularity  and  generates  a  regular  integral.  The  final  term  in  (55),  J3,  has  a 
form  similar  to  that  of  J2, 


1  1 

J3  =  I  j---du2duy, 

ttl=0M2=«l 


the  Duffy  coordinate  transformation  to  use  is. 


(66) 


Ml  =  TJ(0, 

U2  -  CO. 


(67) 


The  three  separate  Duffy  coordinate  transformations  have  projected  the  relative 
coordinate  integration  in  the  J  integrals  onto  the  same  domain.  This  is  the  first  step  in 
combining  the  three  J  terms  together  so  that  can  be  joined  together  in  one  integral. 


Analytic  Integration  over  ^2 


The  inner  integration  over  the  ^  coordinate  in  (50),  (53)  and  (54)  can  be 
performed  exactly  since  the  integrand  depending  on  these  coordinates  can  be  expressed 
as  a  simple  polynomial.  Expanding  and  summing  the  various  terms  in  the  integrand 
produces  the  following  expression. 
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^  ^5  ^)}  ^15^1  ‘*'^14^2  '*'^13^1^2  ^12^1 '^^11^2  "^^10^1  0^112  +  OfUiU2 

+  «7«I  +^6^2  +^A^2  +^3^2^1+02^2^2  +^1- 

(68) 

The  coefficients  of  the  polynomial  are, 

^15  “  2^1  ^14  =  2^4  a^2  =  2/?2  +  2^3  ^12  =  +  2^7  flji  =  2^^  +  2^g 

^10  ■“  ~P\  Og  =  —^4  flg  =  —^2  ~  ^3  ^7  ~  ^5  ^6  ~  ^6  (^^) 

«5=A  ^4=^3  a3=A  «2=A  a^=pg-Alk^ 

where  the  /?s  are  computed  from  (42).  The  analytic  integration  in  the  three  J  integrals 
require  the  following  six  basic  forms, 


b  #,+rf 


K\a,b,c,d)=  J 

^i=a^2=c 

b  4,+d 

K\a,b,c,d)=  j  \m2d^,= 

^i=a#2=c 

b  (i*d 

KHa,b,c,d)=  I 

b  ^y+d 


b  -a  b-a , ,  . 

- + - id-c) 

2  1  ^ 


i3  3  f  2  2 

o  —a  b  —a  .  , 
—  (r/-c) 

b^-a^  b^-a^ 


id-c) 


|,=a42=c 

b  |,+rf 


^i=a^2=c: 

*  ^1+^  l4 

1 

^,=n^2=c 


K  {a,b,c,d)=  \  \^ld^2d^,=^ - - f-  +—(^b-a) 


K\a,b,c,d)=  j  ^^^^^d^^d^,=^--f-  +  dK\aAc,d)  +  ^^^K\aAc,d){^^^ 


(70) 


(71) 


(72) 


(73) 


(74) 
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The  inner  integration  in  the  Ji  integral  is, 

'-"i 


j  J{*(f,5)  +  *«+5-5)te4f, 


(76) 


^1=0  #2=0 


let  this  be  denoted  as  L,(w,,M2)and  because  it  can  be  integrated  analytically,  be  re-written 
as. 
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(0,1  -  w,  ,0,0)  +  (0,1  -  w,  ,0,0)  +  (0,1  -  m,  ,0,0)  + 

(^12  +^^5^1  +aJ^U2)K^{^^--u^S>,0)  +  {a^^+a^u^+a2U2)K^{Q,\-u^,0,Q)  + 

(a,owf  +  a^ul  +  tTgi/jWj  +  +  tZgWj  +  a^  )A!’®(0,1  -  ,0,0). 


Likewise,  denote  the  J2  inner  integration  as  L2(Mi,M2)  which  is  equal  to. 


(77) 


a^^K'^{-u^,\  -  w„-W2,0)  +  a^^K\-u^,\  -  w,  ,-^2,0)  +  0,3^^  (-^2,1  -  u^-u^,^)  + 
(«i2  +  «5“i  ^a^u^)K\-u^,\-u^-u^,Qi)-{-{ai^+a^u^  +a2U2)K^{-u^,\-u^-u^,Q)+ 
(a,QMf  +agMji/2  +a.jU^+a^U2  +aj)jK’°(— 1/2,!- Wi,-M25^X 


and  the  J3  inner  integration  as  Zg(Wi,M2)  which  is  equal  to. 


(78) 


(“2  -  “l  4  -  «i  ,0,  M2  -  Wi )  +  (“2  -  “l  4  -  ,0,  Mj  -  M, )  +  Mi3^*("2  “  ,0,  M2  -  M, ) 

+  (Mj2  +  MjM,  +  a^Mj  )Ar‘  (M2  -  M,  ,1  -  Mj  ,0,  M2  -  M, )  +  (Mj  1  +  M3M,  +  £22 «2  )^^(“2  “  «i  4  “  Wj  ,0,  M2  -  M, ) 

+  (MioM^  +  MgMj  +  M8MjM2  +  a7Mi  +  MgMj  +  Mj  )Ar®  (M2  “  M,  ,1  -  Mj  ,0,  Mj  -  M, ). 

(79) 


This  allows  the  original  singular  integral  to  written  in  compact  form  as  a  regular  integral 
of  the  following  form. 


\  \  th 

inf 

^1=0^2=<^7i=0'72^ 


R(^,n) 


1  ^  ^ 

t|  J  {L^><g{(0,(0r])+L2>^g{a}ri,(0(T]-\))+L2^g{<0T],(0)}(0dr]d(o. 

<y=0  7=0 


(80) 


In  (80)  the  appropriate  Duffy  transformations  are  used  in  each  part  of  the  integrand  to 
map  the  relative  coordinate  m  =(mi,M2)  to  {(0,t])  space. 


Analytic  Integration  over  co 


The  integrand  in  (80)  contains  terms  of  the  form,  and  can  be 

integrated  analytically  with  respect  to  co.  The  algebraic  complexity  of  the  analytical 
integration  is  such  that  it  is  prudent  to  utilize  a  modem  computer  algebra  system  (CAS) 
to  perform  this  tedious  work.  In  addition,  these  CAS  tools  can  generate  Fortran  or  C 
code  directly  and  aid  immensely  in  implementing  these  results  into  existing  or  new 
computer  codes. 
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The  remaining  integration  over  rj ,  must  be  performed  numerically.  However,  the 
analysis  presented  in  this  section  has  transformed  the  original  4-dimensional  singular 
integral  into  a  one-dimensional  non-singular  numerical  integration. 

Common  Edge 

The  instance  where  the  two  triangular  facets  share  a  common  edge  is  formulated 
so  the  two  triangles  share  the  vertices,  V,  and  Vj .  The  relative  distance  vector, 


7-r  =  {t], +(%  -^2 +^.)V2 +#,V3 
=  «.(V,  -  -  V3')  +  (V3  -  V3') 


(81) 


is  a  function  of  the  three  coordinates,  m,,«2  and^,.  The  scalar  distance  R  =  \r-r'\,  is 

zero  only  when  Wj  =  Mj  =  #1  =  0-  In  order  to  remove  the  singularity  that  arises  in  this 
case,  a  3  dimensional  Duffy  transform  is  used  in  the  coordinate  system  defined  by 
w,,«2  and^,. 

The  original  4-dimensional  integral  (9)  in  simplex  coordinates  was  decomposed 
into  a  sum  of  6  separate  integrals. 


^1=0^2=071=072=0  ^ 


(82) 


In  each  of  the  six  integrals  a  3-dimensional  Duffy  transform  will  be  applied  and  the  six 
separate  integrals  will  be  summed  together  to  complete  the  transformation  of  the  original 
singular  integral  into  a  non  singular  integral  of  lower  order. 

Starting  with  the  integral, 

I  I  I  I  (83) 

«,=-!  #2=-U2 

let  z  =  -u ,  and  noting  that  g(i/,^, )  does  not  depend  on  ^2  >  express  Ej  as, 

E,=  j  I  I  (84) 

Zi=0  22=0^1=Zi  ^2=^2 


The  inner  integral,  over  can  be  performed  analytically  since  h  is  a  simple  polynomial 
in  ^2  •  The  domain  of  integration  of  the  remaining  3  dimensional  integral  is  shown  in 
Figure  4,  it  is  a  tetrahedral  with  one  point  at  the  origin.  Equation  (84)  is  not  in  a  form  for 
a  direct  Duffy  transformation,  and  must  be  altered.  The  order  of  the  integration  can  be 
changed  in  a  cyclic  manner,  with  the  aid  of  Figure  4,  so  that  Ei  can  be  re-written  as. 
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(85) 


Equation  (85)  is  now  in  a  form  for  a  direct  Duffy  transform.  The  required  coordinate 
transform  can  be  written  directly. 


r  ^ 

^1 

= 

^(OX^X^j 

(86) 


The  Jacobian  of  the  transformation  is  Equation  (85)  now  can  be  written  in 

standard  form  where  the  coordinates  of  the  g  and  h  functions  are  expressed  in  terms  of 
the  original  u  and  ^  coordinates, 

111  I 

J  J  J  I  —h{-m^-m^x^,ci},^2)^42^^Xydx2dx^da).  (87) 

O)^  JCj  ==0  X2  =0  ^2 

The  integral  E2  is  treated  in  a  similar  fashion, 

^2=  J  1  J  \\K4.'u)g{^A)d^2d^xdu2du^,  (88) 

tt2==0  ^2=0 

where  the  substitution,  r  =  +  Mj  ,  produces. 


The  order  of  integration  is  altered  similar  to  (85), 

^2=1  J  I  g(«>0  I 

r=OMi=O«2==0  52=^ 

The  3-dimensional  Duffy  coordinate  transformation  for  this  integral  is. 


r  ru  ^ 

Ml 

= 

cox^ 

^(OX^  X2^ 

(89) 


(90) 


(91) 


Using  (91)  the  final  form  of  E2  is. 
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1  1  1 


1  1  1  I 


fl;=0jci  =0x2=0 


fj=o 


The  third  integral, 

0  l+«i  I  ^i-(tt2--Wl)  ^  __ 

ttl=-l  ^2"^ 

is  transformed  by  substituting,  z,  =  -u^,Z2  =  ^2  > 


(92) 

(93) 


1  l~2j  1  ^j-(22+rj)  - 

^3=  j  J  1  I  -^K4,z)d^2d4^dz2dz^,  (94) 

2i  =0  22  =0  =22  +2j  ^2  ^ 


and  r  =  Zj  +  Z2 , 


1  1  1  iHh+h) , 

^3=1  I  j  J  -Kzur,<^lA)^^2d^xdrdz^.  (95) 

z,=0r=zi^,=r  ^2=0 

The  integration  domain  for  (95)  is  shown  in  Figure  5,  and  it  is  slightly  different  than  the 
domain  of  Ei.  The  domain  is  a  tetrahedral  with  one  comer  at  the  origin  but  it  is  oriented 
along  the  r  axis.  The  cyclic  reordering  of  the  integration  over  this  domain  lets  E3  be  re¬ 
written  as. 


I  4,  r 

£3=  j  J  j  g(^i.^^i)  J  -7K^ur,^i,42)d42dzidrd^^. 

^i=0r=0z,=0  ^2=0  ^ 

The  3-dimensional  Duffy  coordinate  transformation  for  this  integral  is. 


r 

COX^ 

^0)X^  X2J 

(96) 


(97) 


E3  in  final  form,  written  with  g  and  h  as  functions  of  the  original  simplex  and  relative 
coordinates,  becomes, 

J  j  j  g{-(m^X2,co{\- x^),co)  J  —h{-m^X2,(o{\- x^),(o,^2)d^2^^^\d^2dx\d^ '  (98) 

^-Oxi  =0X2=0  ^2=0 


The  remaining  integrals  are  treated  in  a  similar  fashion  and  the  results  are  summarized 
here. 


Evaluation  of  Singular  Electric  Field  Integral  Equation  matrix  elements 


19 


1  1  1 

-fax,  ,-fax,  jCj  (fa,  jc,  ,  )fa^x,f&2^i<^^ ,  (99) 

<w=0  Jifi=0  JC2=0 

1  I  1 

=  J  j  ^g{m]^,m^x2,a){\-x^))h^{m,x^,x2)ap‘x^dx2clx^d(0 ,  (100) 

<»=0  xi=0  JC2=0 

I  1  1 

=  I  J  ^g{-oxCyX2,a){\-x^),coifi^{(o,x^,X2)G)^x^dx2dx^d(o ,  (101) 

d>=0  Jfj=0  J»^=0 

1  1  i 

£’'‘=11  jg(m^X2,cox^{\-X2),a)(l-x^X2)y^‘*((JD,x^,X2)a)^x^dx2clx^do),  (102) 

<ij=0;ci=0  JC2=0 


1  1  1 

£^  =  j  J  jg(-(ax^X2-(ax^,(o)h^(co,Xy,X2)(O^XlCix2dx^do} ,  (103) 

i»=0jci  =0x2=0 

1  1  1 

£®  =  I  J  ^g(o)Xi,axCiX2,6)(l-x^))h^((o,Xi,X2)co^x^dx2dKid(i).  (104) 

d)=0Xi  =0x2=0 


The  common  edge  Galerkin  inner  product  4-dimensional  integral  (82)  can  now  be 
expressed  as  a  single  integral  with  a  regular  integrand. 


I  ^  J'/ 1:  XT 


j  I  J  J  =  J  J  f  '^g,h'  0)^x^dx2dx^d(o .  (105) 


1  1 


^1=0^2=0  71=0172 


d>=0Xi  =0X2=0 


Analytic  Integration  over  ^2 

The  integrand  of  (105)  is  the  product  of  the  Greens  function,  the  Jacobian  of  the 
3-dimensional  Duffy  transform  and  a  series  of  polynomials.  The  exact  form  of  the 
polynomials,  h‘ ,  are  derived  from  the  following  analytic  integrals. 

fy(l-Xj+XiX2)  ^ 

h^{(o,x^,X2)=  I  — /?(-fux,,-fax,jC2,fU,^2)^^2  (106) 

c  4 

h^{0),X^,X2)=  j  -//(fUX„fUX,.X2,fy(l-JC,),^2¥4  (107) 

#2=0  ^ 

fW(l-Xi)  ^ 

h\(o,x^,X2)=  j  -/z(-fax,X2,fy(l-jr,),fO,4)f/4 

#2=0  ^ 


(108) 
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6){\~X^X2)  j 

h\co,x^,x^)=  j  -h{cox^x^,cox^{\-x^\o){\-x^X2),^2¥^^ 


(109) 


42^(OXi{\-X2) 


h\(0, X^,X2)=  J  , OKiX^ , cy(l  -  or, ), 4 )c/4 

^1=0  ^ 


The  integrand  of  the  h'  expressions  is  given  by. 


(110) 


(111) 


/z(Mj,  Pj’Pk 

+  +  AI2  +  A  («1  +  ^1 )  +  A  ("2  +  4  )  +  A  -  4  / 


In  general,  the  h*  can  be  written  in  this  form, 

A'  (62,  X, ,  Xj )  =  J^(w„  W2 , 4 , 4  )d^2 

a 

h 

~  J®i4  "^^24  "^^3^4 


(113) 


=  {b-d) 


6  +  a 

"T” 


where. 


«i  =  A 

®2  =  A“i  +  A“2  +  (  a + a  )^i  +  A + A 

^3  =  ( A^i  +  A  )“i  +  ( A^i + A )“2  +  ( A + A  )4 + A'^i^  +  A  -  4  / 


(114) 


the  y^^are  given  in  (42).  In  (105)  the  integration  over  6)  can  be  performed  analytically 

because  the  integrand  contains  terms  of  the  form,  .  The  final  result  is  that  2  of 

the  integrations  in  the  original  4  dimensional  integral  can  be  performed  analytically. 


Common  Vertex 

The  case  where  the  two  triangles  share  a  common  vertex  is  a  less  complex 
problem.  If  the  common  vertex  is  chosen  to  be  Vi  then  this  point  is  mapped  to  the  origin 
in  the  simplex  coordinate  space.  The  singularity  arising  from  the  Greens  function  spatial 
distance  term. 
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expressed  in  the  simplex  coordinates  becomes, 

f,  )V,  +«,  - 1,  )V, + -  (n,  - 1?,  )v;  -  V,  V,'.  ( U5) 

Equation  (1 15)  is  zero  only  for. 


Vi  =72  =  ^i  =  4  =0. 


which  is  the  origin  of  the  4  dimensional  simplex  coordinate  space.  The  integral. 


i  1 1 1 

6=0&=0 ia=o  1/2=0  ; 


(116) 


with  the  isolated  singularity  at  the  origin  can  be  regularized  directly  by  using  a  4 
dimensional  generalized  Duffy  transformation  to  eliminate  this  singularity.  In  order  to 
accomplish  this,  the  order  of  the  integration  in  (1 16)  is  first  changed  slightly. 


15  1'?! 


i  ( I  =  ^  { !  f 


^J=:0  7/1=0  ^2=0  72=^ 


/(47) 


(117) 


and  the  domain  split  into  two  parts, 

l  1  4  Vi 


n  n  = 


#1=0  7/i=0  #2=0  7/2=0 

MM 


1  1  4  7i 


1  1  1  ]  i  1  ] 

f,=o,a=o #2=0 1/2=0  ’7  J  Si-Om=?i  ?2=0'?2=o  ; 


(118) 


In  the  second  integral. 


nil  ^^^>iidi,dn,d4,  (119) 

#,=Ov.=fifz=o,ft=o'^V^’7J 


reverse  the  order  of  the  77,,^,,  integration. 


111!  lM^v.dkdn,d4,  =  I  1  n  j§^r,,d4,d4dr,,.  (120) 


#1=0  7/l=#l  #2=0  7/2=0 


7l=0#,=0  #2  =0  7/2=0 


which  is  identical  to  the  first  integral  under  a  coordinate  interchange,  (^  ,;/)=:>  (77  ,^ ). 
The  common  vertex  integral  can  now  be  written  as. 
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i  I  [  I  fin 


1  ^1 


^j=0^2=0  7i=0%=0  ‘ 


4^1=0  r/,^  ^2=0  %=0 


(121) 


The  Duffy  transformation  needed  to  eliminate  the  singularity  at  the  origin  in  this  integral 
has  the  form. 


0)  ^ 

^  COZ2  ^ 

^'52  > 

9 

K^2j 

^CO  Z2Z-^j 

(122) 


This  produces  the  following  result. 


1 1  f  I  n  f  n 

^,=o^2=o^=o%=o 


mm. 


1111 


fi)=0Zj  =0  22=0^3=0 


m.n) ,  /(7.{) 

Rm) 


dz^dz2dZydco. 


(123) 


The  integrand  has  been  left  written  in  terms  of  the  simplex  coordinates  because  of  the 
coordinate  interchange.  To  simplify  the  integrand  is  expressed  as. 


m_m 

mm 


m-ptm 

4 


gi.^m=^mmg{^m- 


(124) 


The  first  term,  h{^,T]),  can  be  expressed  as  a  polynomial, 

mm = midi  ^m<ni +n#2^i +^4724 +«5^i  +«6^2 +^7^1 +a8?72 +^^9  (ns) 

The  Greens  function  depends  only  on  the  distance, 

r-f'  =  (7,  -^.)V,  +te  +^,V,-(7,  -7,)V'  -7,v;, 

(126) 

=  5#i+*4  +  c;7, +r/72. 


where. 


5  =  V2-V,+V3,6  =  -V2,c  =  V,-V'-V;,  J  = 
The  scalar  distance  function  is  expanded  in  a  polynomial, 

mm  =  ^|{d^l  +  ^^2  +  2  Vi  +  ^^2  )■  (5^1  +  ^4  +  c  7i  +  drj2 ), 


(127) 


(128) 


in  compact  form. 
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2  2 


V/2 


The  form  with  the  coordinate  exchange  is. 


EE 

VM  >1 


(129) 


2  2 


vl/2 


II  +  PiJjVi  + 

Vi=i  >1 


(130) 


The  Duffy  coordinate  transform  introduces  a  6>  into  each  of  the  simplex 
coordinates.  The  distance  function  polynomial  in  the  Greens  function  allows  a  single  co 
to  be  brought  out  from  under  the  square  root.  The  Jacobian  of  the  transform  cancels  this 
term  from  the  denominator  of  the  Greens  function  and  eliminates  the  singularity. 

The  original  integral  over  the  pair  of  triangular  facets  with  a  common  vertex  is 
now  written  as. 


tJ  J  J  \^i^^V)gi^,Jl)+Kn4)sinA))G>^Z2dzjdz^dz^d(o.  (131) 


The  integration  over  (o  can  be  performed  analytically  because  the  integrand  contains 

terms  of  the  form,  The  remaining  integrations  must  be  performed 

numerically. 

Summary  and  Conclusion 

The  Galerkin  solution  to  the  electric  field  integral  equation  using  RWG  basis 
functions  has  and  will  continue  to  be  popular  for  simulating  electromagnetic  phenomena. 
Fundamental  to  achieving  an  accurate  result  is  the  numerical  accuracy  of  the  impedance 
matrix  elements  that  are  used  in  computing  the  solution.  The  material  presented  here 
shows  how  to  improve  the  computation  of  matrix  elements  that  involve  singular  surface 
integrals  by  removing  the  singularity  in  the  integrand  through  a  series  of  coordinate 
transforms.  The  formulation  provided  an  additional  benefit  by  reducing  the  dimension  of 
the  4-dimensional  integrations  that  result  from  the  Galerkin  approach  by  1,  2  or  3 
depending  on  the  degree  of  overlap  between  the  two  triangular  patches. 

This  analysis  can  be  extended  to  curved  geometries,  quadrilateral  patch  geometry 
and  other  integral  equation  operators  encountered  in  electromagnetics.  The  extension  of 
these  methods  to  hypersingular  integral  equation  operators  will  provide  a  complete 
foundation  for  highly  accurate  method  of  moments  numerical  solutions. 

The  material  presented  here  is  a  modernization  of  the  classic  RWG  method.  The 
next  step  in  this  modernization  process  is  to  utilize  computer  algebra  software,  e.g.  Maple 
or  Mathematica,  to  complete  the  tedious  algebraic  manipulations  and  analytic 
integrations.  The  ultimate  goal  is  to  generate  modular  Fortran90/95  code  that  could  be 
used  as  a  foundation  for  modem  method  of  moments  (MoM)  electromagnetic  analysis 
tools. 
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Figure  1.  Rao,  Wilton,  Glisson  (RWG)  vector  basis  function.  The  basis  function 
consists  of  a  pair  of  triangular  patches,  with  areas  Ap  and  Aq,  which  share  a  common  edge 
with  length  /„.  The  expansion  of  the  unknown  surface  electric  current  in  terms  of  the 
vector  function  /(r )  is  associated  with  the  n*  edge. 
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Figure  2.  The  domain  of  integration  of  the  4  dimensional  Galerkin  irmer  product  integral 
expressed  in  terms  of  relative  and  simplex  coordinates.  This  is  used  for  interchanging  the 
order  of  integration  in  the  integrals  l'  through 


Figure  3.  The  domain  of  integrati 
and  J3,  that  result  from  the  commoi 
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Zl 


Figure  4.  The  domain  of  integration  of  the  integral  Ei  for  the  common  edge  case.  A  unit 
cube  is  shown  for  reference,  the  red  lines  define  a  tetrahedral  region  which  border  the 
domain  of  integration  of  the  Ei  integral. 
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Figure  5.  The  domain  of  integration  of  the  integral  E3  for  the  common  edge  case.  A  unit 
cube  is  shown  for  reference,  the  red  lines  define  a  tetrahedral  region  which  border  the 
volume  of  integration  of  the  E3  integral. 


